Approximation of conformal mappings by circle 

patterns 

Ulrike Bucking * 

Abstract 

A circle pattern is a configuration of circles in the plane whose combi- 
natorics is given by a planar graph G such that to each vertex of G there 
corresponds a circle. If two vertices are connected by an edge in G then 
the corresponding circles intersect with an intersection angle in (0, tt) and 
these intersection points can be associated to the dual graph G* . 

Two sequences of circle patterns are employed to approximate a given 
conformal map g and its first derivative. For the domain of g we use 
embedded circle patterns where all circles have the same radius e n > 
for a sequence e„ — > and where the intersection angles are uniformly 
bounded. The image circle patterns have the same combinatorics and 
intersection angles and are determined from boundary conditions (radii 
or angles) according to the values of g' (\g'\ or arg g'). The error is of order 
l/y/— loge n . For quasicrystallic circle patterns the convergence result is 
strengthened to C°°-convergence on compact subsets and an error of order 



1 Introduction 

Conformal mappings constitute an important class in the field of complex analy- 
sis. They may be characterized by the fact that infinitesimal circles are mapped 
to infinitesimal circles. Suitable discrete analogs are of actual interest in the 
area of discrete differential geometry and its applications, see [BSSZ08j . 

Bill Thurston first introduced in his talk |Thu85j the idea to use finite circles, 
in particular circle packings, to define a discrete conformal mapping. Remember 
that an embedded planar circle packing is a configuration of closed disks with 
disjoint interiors in the plane C. Connecting the centers of touching disks by 
straight lines yields the tangency graph. Let and ^2 be two circle packings 
whose tangency graphs are combinatorially the same. Then there is a mapping 
g<# : ^1 — * ^2 which maps the centers of circles of ^\ to the corresponding- 
centers of circles of %?2 and is an affine map on each triangular region corre- 
sponding to three mutually tangent circles. Various connections between circle 
packings and classical complex analysis have already been studied. A beautiful 
introduction and surway is presented by Stephenson in Stc05 . In particular, 
several results concerning convergence, i.e. quantitative approximation of con- 
formal mappings by g v , have been obtained, see [RS871 IClT92l IHS961 IHS98] . 
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The class of circle patterns generalizes circle packings as for each circle pack- 
ing there is an associated orthogonal circle pattern. Simply add a circle for each 
triangular face which passes through the three touching points. To define a cir- 
cle pattern we use a planar graph as combinatorial data. The circles correspond 
to vertices and the edges specify which circles should intersect. The intersection 
angles arc given using a labelling on the edges. Thus an edge corresponds to 
a kite of two intersecting circles as in Figure [1] (right). Moreover, for interior 
vertices the kites corresponding to the incident edges have disjoint interiors and 
their union is homcomorphic to a closed disk. Similarly as for circle packings, 
there are results on existence, rigidity, and construction of special circle pat- 
tern. See for example |Riv94( ISch97j IAB00) IBSM) IBH03j . where some of the 
references use a generalized notion of circle patterns. 

Given two circle patterns ^\ and ^2 with the same combinatorics and inter- 
section angles, define a mapping g^ : ^ — > similarly as for circle packings. 
Namely, take g<g to map the centers of circles and the intersection points of ^1 
corresponding to vertices and faces of G to the corresponding centers of circles 
and intersection points of ^2 and extend it to an affine map on each kite. 

For a given conformal map g we use an analytic approach and specify suit- 
able boundary values for the radius or the angle function according to \g'\ or to 
arg</ respectively in order to define the (approximating) mappings g^g. Gener- 
alizing ideas of Schramm's convergence proof in |Sch97j we obtain convergence 
in C 1 on compact sets if we take for a sequence of isoradial circle patterns 
(i.e. all radii are equal) with decreasing radii e n — > which approximate the 
domain of g. Furthermore, we assume that the intersection angles are uniformly 
bounded away from and tt. Note in particular that the combinatorics of the 
circle pattern < ^ n ' may be irregular or change within the sequence. Thus our 
convergence results applies to a considerably broader class of circle patterns as 
the known results of Schramm [Sch97| , Matthcs [Mat05| , or Lan and Dai [LD07] 
for orthogonal circle patterns with square grid combinatorics. 

The main idea of the proof is to consider a "nonlinear discrete Laplace 
equation" for the radius function. This equation turns out to be a (good) 
approximation of a known linear Laplace equation and can be used in the case 
of isoradial circle patterns to compare discrete and smooth solutions of the 
corresponding elliptic problems, that is the logarithm of the radii of ^2"' and 
log \g'\ — Re(log g'). Also, we obtain an a priori estimation of the approximation 
error of order 1 / y/\og{l/e n ). 

If the circle patterns additionally have only a uniformly bounded num- 
ber of different edge directions, then the corresponding kite patterns are qua- 
sicrystallic rhombic embeddings and the circle patterns ^ are called qua- 
sicrystallic BMS05 . For such embeddings we generalize an asymptotic devel- 
opment given by Kenyon in Kc n02j of a discrete Green's function. Also, using 
similar ideas as Duffin in |Duf53j , we generalize theorems of discrete potential 
theory concerning the regularity of solutions of a discrete Laplace equation. We 
then use these results together with the C 1 -convergence for isoradial circle pat- 
terns and prove C°°-convergencc on compact sets for a class of quasicrystallic 
circle patterns. In this case the approximation error is of order e n (or e\ for 
square grid or hexagonal combinatorics and regular intersection angles). The 
proof generalizes a method used by He and Schramm in |HS98j . 

The article is organized as follows: First we introduce and remind the ter- 
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urinology and some results on circle patterns in Section^ focusing in particular 
on the radius and the angle function. In Section [3] we formulate and prove the 
theorems on C^-convergence for isoradial circle patterns. After a brief review 
on quasicrystallic circle patterns in Section 0] we state and prove in Section [5] 
the theorem on C°°-convergence. The necessary results on discrete potential 
theory arc presented in Appendix [A"l An extended and more details version of 
the results can be found in |Biic07j . 

The author would like to thank Alexander I. Bobenko, Boris Springborn, 
and Yuri Suris for various discussions and helpful advice. 

2 Circle patterns 

To define circle patterns we use combinatorial data and intersection angles. 

The combinatorics are specified by a b-quad-graph 3>, that is a strongly 
regular cell decomposition of a domain in C possibly with boundary such that 
all 2-cells (faces) are embedded and counterclockwise oriented. Furthermore all 
faces of £F are quadrilaterals, that is there are exactly four edges incident to 
each face, and the 1-skeleton of ^ is a bipartite graph. We always assume that 
the vertices of & are colored white and black. To these two sets of vertices we 
associate two planar graphs G and G* as follows. The vertices V(G) are all 
white vertices of V(£)). The edges E(G) correspond to faces of that is two 
vertices of G are connected by an edge if and only if they are incident to the 
same face of 3). The dual graph G* is constructed analogously by taking for 
V(G*) all black vertices of V(S#). £F is called simply connected if it is the cell 
decomposition of a simply connected domain of C and if every closed chain of 
faces is null homotopic in 3. 

For the intersection angles, we use a labelling a : F(2>) — > (0, 7r) of the faces 
of S>. By abuse of notation, a can also be understood as a function defined 
on E(G) or on E(G*). The labelling a is called admissible if it satisfies the 
following condition at all interior black vertices v £ Vi n t(G*): 



Figure 1: Left: An example of a b-quad-graph 3) (black edges and bicolored 
vertices) and its associated graph G (dashed edges and white vertices). Right: 
The exterior intersection angle a of two intersecting circles and the associated 
kite built from centers and intersection points. 





(1) 
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(a) A regular square (b) A regular hexagonal (c) A rhombic embedding (a part of a 
grid circle pattern. circle pattern. Penrose tiling) and a corresponding 

isoradial circle pattern. 



Figure 2: Examples of isoradial circle patterns. 

Definition 2.1. Let ^ be a b-quad-graph and let a : E(G) — > (0,7r) be an 
admissible labelling. An (immersed planar) circle pattern for @ (or G) and a 
are an indexed collection c (o = {C z : z G V(G)} of circles in C and an indexed 
collection = {K e : e G E(G)} = {K f : f G F{@)} of closed kites, which all 
carry the same orientation, such that the following conditions hold. 

(1) If z\,z-2 G V(G) are incident vertices in G, the corresponding circles C Zl , C Z2 
intersect with exterior intersection angle a([2ii, 22])- Furthermore, the kite 
-fTr Zlj22 ] is bounded by the centers of the circles C Zl , C Z2 , the two intersection 
points, and the corresponding edges, as in Figure[T] (right). The intersection 
points are associated to black vertices of V{S#) or to vertices of V(G*). 

(2) If two faces are incident in then the corresponding kites have one edge 
in common. 

(3) Let /i,...,/„ G F{&) be the faces incident to an interior vertex v G 
Vint(^). Then the kites K / t , . . . , Kf n have mutually disjoint interiors. 
Their union Kf 1 U • • • U Kf n is homeomorphic to a closed disc and con- 
tains the point p(v) corresponding to v in its interior. 

The circle pattern is called embedded if all kites of Jtf have mutually disjoint 
interiors. It is called isoradial if all circles of ^ have the same radius. 

Note that we associate a circle pattern ^ to an immersion of the kite pattern 

corresponding to & where the edges incident to the same white vertex are 
of equal length. The kites can also be reconstructed from the set of circles 
using the combinatorics of G. Note further, that there are in general additional 
intersection points of circles which are not associated to black vertices of V{&). 
Some examples are shown in Figure [H 

There are also other definitions for circle patterns, for example associated 
to a Delaunay decomposition of a domain in C. This is a cell decomposition 
such that the boundary of each face is a polygon with straight edges which is 
inscribed in a circular disk, and these disks have no vertices in their interior. The 
corresponding circle pattern can be associated to the graph G* . The Poincare- 
dual decomposition of a Delaunay decomposition with the centers of the circles 
as vertices and straight edges is a Dirichlet decomposition (or Voronoi diagram) 
and corresponds to the graph G. 

Furthermore the definition of circle patterns can be extended allowing cone- 
like singularities in the vertices; see [BS04] and the references therein. 
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2.1 The radius function 

Our study of a planar circle pattern ^ is based on characterizations and prop- 
erties of its radius function r^ = r which assigns to every vertex z <E V(G) 
the radius r<g(z) — r(z) of the corresponding circle C z . The index 'rf will be 
dropped whenever there is no confusion likely. 

The following proposition specifies a necessary and sufficient condition for a 
radius function to originate from a planar circle pattern, see |BS04] for a proof. 
For the special case of orthogonal circle patterns with the combinatorics of the 
square grid, there are also other characterizations, see for example |Sch97j . 

Proposition 2.2. Let G be a graph constructed from a b- quad- graph £F and let 
a be an admissible labelling. 

Suppose that ^ is a planar circle pattern for & and a with radius function 
r = TV . Then for every interior vertex Zq € Vint (G) we have 



and the branch of the logarithm is chosen such that < fe{x) < n. 

Conversely, suppose that @ is simply connected and that r : V(G) — > (0, oo) 
satisfies @ for every z £ Vint{G). Then there is a planar circle pattern for 
G and a whose radius function coincides with r. This pattern is unique up to 
isometries of C. 

Note that < 2,f a ([z,z }) Qogr(z)— logr(zo)) is the angle at zo of the kite with edge 
lengths r(z) and r(z ) and angle a([z, z ]), as in Figure [T] (right). Equation 
is the closing condition for the chain of kites corresponding to the edges incident 
to zq which is condition (3) of Definition 12.11 

For further use we mention some properties of fg, see for example |Spr03| . 



Lemma 2.3. (1) The derivative of fg is fg{x) = 2(cosh'"-cos W) > ^' ^° f 9 * s 



strictly increasing. 

(2) The function fg satisfies the functional equation fg{x) + fe(~x) = n — 9. 

(3) For < y < ir - 9 the inverse function of fg is fg 1 (y) = log sil ^)| e) ■ 

Remark 2.4. Equation (0) can be interpreted as a nonlinear Laplace equation 
for the radius function and is related to a linear discrete Laplacian which is 
common in the linear theory of discrete holomorphic functions; see for exam- 
ple DufGS, McrOl, BMS05 for more details. This can be seen as follows. 

Let G be a planar graph and let a be an admissible labelling. Assume there 
is a smooth one parameter family of planar circle patterns ^ e for G and a with 
radius function r e for e £ (—1,1). Then for every interior vertex zq with incident 
vertices Z\, . . . ,z m and all e G (—1, 1) Proposition 1 2 . 21 implies that 




where 



fe{x) ■= — log 



2 _ e x+i9 ' 



m 



^ 2 /«([^o])( lo g r £( z j) - logr £ (z )) = 2tt. 
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Differentiating this equation with respect to e at e = 0, we obtain 



77} 



^2 2 fL([z j ,z ])( l °Sro( z j) - logr (^o))(f(^) - v(z )) = 0, 



(3) 



3=1 



where v(z) = logr £ (2;)| £= o- Thus v satisfies a linear discrete Laplace equation 
with positive weights. Lemma 12.31 and a simple calculation show that 



Here z\, Zi £ V(G) are two incident vertices which correspond to the centers of 
circles p(z\),p(z2) of the circle pattern ^o- The two other corner points of the 
same kite are denoted by p(v{),p(v2)- 

In analogy to smooth harmonic functions, the radius function of a planar 
circle pattern satisfies a maximum principle and a Dirichlct principle. 

Lemma 2.5 (Maximum Principle). Let G be a finite graph associated to a b- 
quad-graph as above with some admissible labelling a. Suppose %f and 'ta* are 
two planar circle patterns for G and a with radius functions r%>,r<g» : V{G) — > 
(0, oo). Then the maximum and minimum of the quotient rv /r^* is attained at 
the boundary. 

A proof can be found in Hc99l Lemma 2.1]. If there exists an isoradial 
planar circle pattern for G and a, the usual maximum principle for the radius 
function follows by taking r^g- = 1. 

Theorem 2.6 (Dirichlct Principle). Let Si be a finite simply connected b-quad- 
graph with associated graph G and let a be an admissible labelling. 

Let r : Vg(G) — * (0, oo) be some positive function on the boundary vertices 
of G. Then r can be extended to V(G) in such a way that equation ^ holds 
at every interior vertex z £ Vi n t(G) if and only if there exists any circle pattern 
for G and a. If it exists, the extension is unique. 

By lack of a good reference we include a proof. 

Proof. The only if part follows directly from the second part of Proposition ^. 21 
To show the if part, assume that there exists a circle pattern for G and a 
with radius function R : V(G) — > (0, oo). A function k : V(G) — > (0, oo) which 
satisfies the inequality 



at every interior vertex z £ Vi n t{G) will be called subharmonic in G. Let b be 
the minimum of the quotient r/R on Vq(G) and let k,\ be equal to r on Vg(G) 
and to bR on Vi n t{G). Then k\ is clearly subharmonic. The maximum of K\/R 
is attained at the boundary, which is a simple generalization of the Maximum 
Principle 12.51 Let r* be the supremum of all subharmonic functions on G that 
coincide with r on Vq{G). Thus r* is bounded from above by the maximum of 
r/R on Vg(G) which is finite. One easily checks that r* satisfies condition 
The uniqueness claim follows directly from the Maximum Principle 12.51 □ 



2 f' a {[z u z 2 ]){^ r o{ z i) - logr (z 2 )) 



P(vi)-p(v2) 
p(zi) -p(z 2 ) 



(4) 
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2.2 The angle function and relations to the radius function 

Similarly as for polar coordinates of the complex plane, a suitably defined angle 
function can be interpreted as a "dual" to the radius function. We focus on con- 
nections between these functions and on a characterization for angle functions 
of circle patterns similar to Proposition ^. 21 

Let G be a finite graph associated to a b-quad-graph S and let a be an 
admissible labelling. Denote by E(3>) the set of oriented edges, where each 
edge of E(3i) is replaced by two oriented edges of opposite orientation. Let 
be a planar circle pattern for S and a. Define an angle function ip^ — ip on 
E(S) as follows. Denote by p(w) the point of ^ corresponding to the vertex 
w E V(S). For e = zv £ E(£)) set ip(e) = arg(p(w) — p(z)) to be the argument 
of p(v) — p(z), that is the angle between the positively oriented real axis and the 
vector p(v) — p(z). As this argument is only unique up to addition of multiples 
of 2tt, we will mostly consider ip G R/(27rZ). But note that e lv ^ is well defined. 
Furthermore 

(p(e) — <p(—e) = 7r (mod 2-k). (6) 

Our choice of the angle function leads to the following connections between 
radius function r and angle function ip for a planar circle pattern for and a. 

Let / G F{&) be a face of &. Without loss of generality, we assume that 
the notation for the vertices and edges of / is taken from Figure [3] (left). More 
precisely, the white vertices z_ and z + of / and the black vertices u_ and v + are 
labelled such that the points z_ , V— , z + , v + appear in this cyclical order using 
the counterclockwise orientation of /. Furthermore the edges are labelled such 
that e*i = Z-V-l, e*2 = Z-vl, — e*3 = i>_z+, — e*4 = v+z+. Then the following 
equations hold. 

p(ei) - </>(-e 3 ) = a(f) - w + 2f a{f) (log r(z + ) -logr(z_)) (mod 2tt) (7) 
ip(-S 4 ) - <p(S 2 ) = a{f) -7T + 2/„ (/) (log r(z + ) -logr(z_)) (mod 2tt) (8) 
(p(ei) - (p(e 2 ) = 2/ a(/) (logr(z + ) - logr(z_)) (mod 2tt) (9) 

<p(-e 4 )-<p(-e 3 ) = -2f a(f) (logr(z-)-logr(z + )) (mod 2^) (10) 

<f(-e3) - <f(e2) = 7r - a{f) (mod 2tt) (11) 

tp(-ei) - v?(ei) =a(/)-7r (mod 2tt) (12) 

Additionally, the angle function <p satisfies the following 

Lemma 2.7 (Monotonicity condition). Let z G ^(G) 6e a w/iiie vertex and 
let ei, . . . , e„ &e t/ie sequence of all incident edges in which are cyclically 

ordered respecting the counterclockwise orientation of the circle C z . Then the 
values of <p G R can be changed by suitably adding multiples of 2tt such that 
p>(ej) is an increasing function of the index j and if z is an interior vertex, then 
also ip(ej) — f(ei) < 2ir for all j = 1, 



, n. 



The following theorem is useful to compare two circle patterns with the same 
combinatorics and intersection angles. 

Theorem 2.8. Let & be a b-quad-graph and let a be an admissible labelling. 

Let c ta and be two planar circle patterns for &! and a with radius functions 
fcg = r and r^ = f and angle functions tp<g = if and ip^, = (p respectively. 
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Figure 3: Left: A face of & with oriented edges. Right: An interior intersection 
point with its neighboring faces. 



Then the difference tp — tp gives rise to a function 8 : V(G*) —> R such that the 
following condition holds on every face f € F(S>). 

0°g + io § - 2/-w 0°g = - *(«-) 

(13) 

_ffere we have defined w : V(G) — > R, u>(z) = f(z)/r(z) and the notation is 
taken from Figured (left) as above. 

Conversely, assume that ^ is a planar circle pattern for Si and a with radius 
functions r^g = r and angle function tp<g = tp. Let 8 : V(G*) — > R and w : 
V(G) —* R+ be two functions which satisfy equation Q13p for every face f € 
F{$>). Then (rw) and (ip + 8) are the radius and angle function of a planar 
circle pattern for £1 and a. This circle pattern is unique up to translation. 

Proof. If and ^ are two planar circle patterns for @ and a, equations ((6]), 
(fTTj) , and (fT2|) imply that the difference tp — ip is constant for all edges incident 
to any fixed black vertex v £ V(G*). Therefore 8 (mod 2n) is well defined on 
the intersection points and encodes the relative rotation of the star of edges at 
v. Also equation Q13p holds modulo 2tt. 

To obtain a function 5 with values in R, fix S(vq) G [0, 2tt) for one arbitrary 
vertex i>o € V(G*) (for each connected component of G* ). Define the values of 8 
for all incident vertices in G* by equation (fT3"|) . Continue this construction until 
a value has been assigned to all vertices. This proceedure leads to a well-defined 
function, as by Proposition 12. 21 the sum of the left hand side of equation (fT3"|) is 
zero for simple closed paths in E(G*) around a white vertex. 

To prove the converse claim, observe that if r is a radius function, equa- 
tion (fT5|) implies that (wr) fullfills equation By Proposition 12.21 there is a 
circle pattern with radius function (wr). Adjust the rotational freedom at one 
edge according to (ip + S). Equation (fO|l implies that (tp + 8) is indeed the angle 
function of this circle pattern. □ 

The preceeding theorem motivates the definition of a comparison function 
for two circle patterns with the same combinatorics and intersection angles. 

Let G be a graph associated to a b-quad-graph *2> and let a be an admissible 
labelling. Suppose that ^\ and ^2 are planar circle patterns for $ and a with 
radius functions and r<^ 2 and angle functions tp^ and tp<g 2 respectively. Let 
8 : V(G*) — > R be a function corresponding to tp<g 2 — tp^ as in Lemma 12.81 
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Define a comparison function w : V(&) — ► C by 

w(y) = T^Av)l r ^Av) for y e v(G), 

= e* 5 ^ e S 1 for x G V(G*). 

Note that u>(y) is the scaling factor of the circle corresponding to y £ V(G) 
when changing from the circle pattern c <fi to ^2- w(x) gives the rotation of the 
edge-star at a; € V(G*). Furthermore, w satisfies the following Hirota Equation 
for all faces / £ F{&). 

w(x )w(y )a Q - w(xi)w(y )ai - w{x 1 )w(y 1 )a + w(x )w(yi)a 1 = (15) 

Here xq, x\ £ V(G*) and yo, y\ £ V(G) are the black and white vertices incident 
to / and ao = xq — yo and a\ = x\ — yo are the directed edges. Equation (|15p 
is the closing condition for the kite of ^2 which corresponds to the face /. 

Angle functions associated to planar circle patterns can be characterized in 
a similar way as radius functions are qualified in Proposition ^. 21 

Proposition 2.9. Let S 1 be a b- quad- graph with associated graphs G and G* 
and let a be an admissible labelling. 

Suppose & is a planar circle pattern for & and a with angle function ip = tp<# . 
Then ip satisfies equations ©, Q11JI . (|12|) . the Monotonicity condition \ 2.7\ at 
every white vertex of 3 ' , and the following two conditions. 

(i) Let f be a face of S> and let e\ and e-2. be two edges incident to f and to the 
same white vertex and assume that e\ and e^ are enumerated in clockwise 
order as in Figured (left). Define an angle (3 £ (0,7r) by 

2/3 = (p(e{) - <p(e*!) (mod 2w), 

where the orientation of the edges is chosen such that the vectors point 
from a white vertex to a black vertex, as in Figure^ (left). Then 

/3 + a(/)<7r. (16) 



(ii) For an interior black vertex v £ Vi n t{G*), denote by e\, . . . , e n ,e n +i — e\ 
all incident edges of S> in counterclockwise order and by fj the face of S> 
incident to ej and e-j+i- Denote by e* (j = l,...,n) the edge incident 
to ej and fj which is not incident to v. For j = 1, . . . , n define as above 
j3j £ (0,7r) by 2f3j = ip{€j) — tp(e~j ) (mod 2tt), where we choose the same 
orientation of the edges from white to black vertices as above. Then 



'«(/: 

3=1 



•m(&)=0. (17) 



Conversely, suppose that @> is simply connected and that ip : E(@) — > 
R/(27rZ) satisfies equations ©, ([TTj) , (112p . the Monotonicity condition \2. 7[ 
condition (|16|) at every white vertex of 3$, and condition (|17p at every inte- 
rior black vertex. Then there is a planar circle pattern for 3> and a with angle 
function ip. This pattern is unique up to scaling and translation. 



10 



Ulrike Bucking 



Proof. For a given planar circle pattern equations ([7])- (fT^)) and Lemma |2~T1 hold . 
To show (fT()|) . consider a kite corresponding to a face of S> as in Figure [3] (left). 
Note that /3_ = 2/3 by equation ©, /3„ + /?+ + 2a = 2n, and [3-,f3 + ,a > 0. 
Using notation of Figure[3] (right) we also deduce that f~} (f3j ) = log Tj+\ — log Tj 
for j = I, ... where we identify r\ = r n +i. Now (|17|) follows immediately. 

In order to prove the converse claim, we construct a radius function r : 
V(G) — > (0, oo) and build a circle pattern corresponding to r and p. 

Let z £ V(G) be an interior white vertex. Set r(z) = 1. Consider a neigh- 
boring white vertex z' £ V(G) and the face / £ F(&) incident to z and z' . 
Denote the black vertices of & incident to z and / by v\, V2 such that v\, 
z, V2, z' appear in counterclockwise order along the boundary of /. Define 
ifj- £ (0, 271") by tp- = (p(zvi) — (p(zv2) (mod 2tt). Condition ([TBJ implies that 
r(z') := r(z)exp(/~^(^_/2)) is well defined and positive. We proceed in this 
way until a radius has been assigned to all white vertices. Condition (TT7|) guaran- 
tees that these assignments do not lead to different values when turning around 
a black vertex (see Figure[3] (right) with ip = 2(3). Thus r is uniquely determined 
up to the choice of the initial radius, which corresponds to a global scaling. For 
each face / £ F(3t) construct a kite with lengths r(z + ), r(z_) of the edges inci- 
dent to the white vertices z + , z_ of / respectively and angle a = a(f). Lay out 
one kite fixing the rotational freedom according to ip. Successively add all other 
kites, respecting the combinatorics of By construction and assumptions, at 
every interior vertex the angles of the kites having this vertex in common add 
up to 2tt. Thus we obtain a circle pattern with angle function tp. □ 

3 C 1 -convergence with Dirichlet or Neumann 
boundary conditions 

In this section we state and prove our main results on convergence for isoradial 
circle patterns. We begin with Dirichlet boundary conditions, that is we first 
focus on the radius function with given boundary values. 

Theorem 3.1. Let D C C be a simply connected bounded domain, and let 
W C C be open such that W contains the closure D of D. Let g : W — > C be a 
locally infective holomorphic function. Assume, for convenience, that £ D. 

For n £ IN let 3> n be a b- quad-graph with associated graphs G n and G* and 
let a n be an admissible labelling. We assume that £) n is simply connected and 
that a n is uniformly bounded such that for all n £ IN and all faces f £ F{S> n ) 

\a n (f) - tt/2| < C (18) 

with some constant < C < tt/2 independent of n. 

Let e n £ (0, oo) be a sequence of positive numbers such that e n — > for 
n — ► oo. For each n £ IN, assume that there is an isoradial circle pattern for 
G n and a n , where all circles have the same radius e n . Assume further that all 
centers of circles lie in the domain D and that any point x £ D which is not 
contained in any of the disks bounded by the circles of the pattern has a distance 
less than Ce n to the nearest center of a circle and to the boundary dD, where 
C > is some constant independent of n. Denote by R n = e n and <f) n the radius 
and the angle function of the above circle pattern for G n and a n . By abuse of 
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notation, we do not distinguish between the realization of the circle pattern, that 
is the centers of circles z n , the intersection points v n , and the edges connecting 
corresponding points in 3i n or G n , and the abstract b-quad-graph 3i n and the 
graphs G n and G* . Also, the index n will be dropped from the notation of the 
vertices and the edges. 

Define another radius function on G n as follows. At boundary vertices z G 
V d (G n ) set 

r n {z)=R n {z)\g'{z)\. (19) 

Using Theorem \2.6\ extend r n to a solution of the Dirichlet problem on G n . Let 
zq G V(G n ) be such that the disk bounded by the circle C Zo contains and let 
e = [zq,vq] G E(S> n ) be one of the edges incident to zq such that 4> n (e) G [0, 2ir) 
is minimal. 

Let ip n be the angle function corresponding to r n that satisfies 

tp n (g) = arg(» )) + <£„(e). (20) 

Let ^„ be the planar circle pattern with radius function r n and angle function 
ip n . Suppose that ^ n is normalized by a translation such that 

Pn(vo) = 9(vo), (21) 

where p n (v) denotes the intersection point corresponding to v G V(G*J. For 
z G D set 

g n (z)=p n (w) and q n (z) = ^l e ^M)-^M)) ; 

R n {v) 

where w is a vertex o/ V(G*) closest to z and v is a vertex ofV{G n ) closest to 
z such that [v,w] G E(@ n ). 

Then q n — > g' and g n — > g uniformly on compact subsets in D as n — > oo. 

Remark 3.2. The proof of Theorem 13.11 actually shows the following a priori 
estimations for the approximating functions q n and g n . 

lkn-3'||y(G„)nJ<r < Gi(-log 2 e„)"5 and \\g n -g\\ V (G n )r\K < C 2 {- log 2 e n )~i 

for all compact sets K C D, where the constants Ci,C 2 depend on K, g, D, 
and on the constants of Theorem 13.11 

We begin with an a priori estimation for the quotients of the radius functions. 

Lemma 3.3. For z G V(G n ) set 

h n (z) = log 1^(^)1 , 

t n (z) = log(r n (z) / R n (z)) = log(r„(z)/e„). 

Then 

h n (z)-t n (z) = 0(e n ). 

Here and below the notation Si = 0(s2) means that there is a constant 
G which may depend on W,D,g, but not on n and z, such that \s\\ < Gs 2 
wherever s\ is defined. A direct consequence of Lemma 13.31 is 

r n (z) = R n (z)\g'(z)\+0(e 2 n ). (22) 

Our proof uses ideas of Schramm's proof of the corresponding Lemma in |Sch97j . 
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Proof. Consider the function 

p(z)=t n (z) - h n {z)+(3\z\ 2 , 

where /3 € (0, 1) is some function of e n . We want to choose (3 such that p will 
have no maximum in Vi„t(G n ). 

Suppose that p has a maximum at z G Vi n t(G n ). Denote by z\, . . 
incident vertices of z in G n in counterclockwise order. Then we have 



^n(^j/) ^n(^) — Xj 



i Zm the 
(23) 



for j = 1, 



where 



Xj = Kizj) - h n {z) - (3\ Z] \ 2 + /3|z| 5 



(24) 
0(e n ). 



Since z € Vi n t(G n ), we have \z\ = 0(1) and by assumption z — z 
With (3 e (0, 1) this leads to (3\zj\ 2 - (3\z\ 2 = 0(e n ). Using this estimate and 
the smoothness of Re (log g'), we get Xj = 0(e n ). 

From ([2"3"|) , the definition of t n (z) = \og(r n (z) / R n (z)) = logr„(z) — loge n 
and the monotonicity of the sum in equation (J2|) (see Lemma 12.31 (i)), we get 



^ m 

fa[z,z s )0ogrj^j) - logr„(z)) 

V = 1 =t n {Zi)-t n {z) 



7T < 



^2fa(z,zj) (Xj) 
^=1 



(25) 



Remembering = 0{e n ), we can consider a Taylor expansion of the right hand 
side of inequality (|25[) about in order to make an 0(e^)-analysis. 

Consider the chain of faces fj of 2$ n (j 

1 m) which are incident to z and Zj. The 

enumeration of the vertices Zj (and hence of the 
faces fj) and of the black vertices Vi, . . . , v rn in- 
cident to these faces can be chosen such that fj 
is incident to Wj-i and Vj for j = 1, . . . , m, where 
= v m . Furthermore, using this enumeration 
i{zj — z) and (vj — Vj-i) are parallel, see Figured] 
As each face fj of an isoradial circle pattern is 
a rhombus we can write, using the notation of 




Vj-l 

Figure 4: A rhombic face 
of & n with oriented edges. 



Figure |H 



a-i-i 



and 



Vj-i 



(26) 



Denoting aj — a([z,Zj]), lj 



z\ = 2e„ sin(<x,-/2), and lj = \vj — Vj^i\ 



2e n cos(aj /2), we easily obtain by simple calculations that 

/«. (0) = (tt - a 3 -)/2, 4 (0) = lj/{2lj), fl 3 (0) = 0. 

Taking into account that equation (|2|) holds with R n = e n and using the uniform 
boundedness (|18|) of the labelling a, inequality (j25|) yields 



o< ^4(0)^ + 0(4). 

i=i 



(27) 
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To evaluate this sum, expand 

h n { Zj ) - h n {z) = Re(\ogg'{ Zj ) - logs' (z)) = Re(a( Zj - z) + b( Zj - z) 2 ) + 0{e 3 n ) 
and Xj = h n (Zj) - h n (z) - (3\zj\ 2 + f3\z\ 2 

= Rc(a( Zj -z) + b(zj - z f - 2pz( Zj - z)) - f3l 2 + 0{e 3 n ). 

Noting that f' aj {0)(zj — z) = (vj — Vj-\)/{2i) we get 

m fa — 2j3z m b m \ 

3 = 1 \ 3=1 3 = 1^. ' '/ 



77. 

3=1 

Thus from inequality (|27p . remembering lj = 2e n sin(a :) /2) and lj = 2e n cos((Xj/2), 
we arrive at 

m tci 

< -0e 2 n sin(a 3 -/2) cos(a 3 -/2) + 0(e 3 n ) /Sgsinfo-) < 0(e n ). 

3 = 1 3 = 1 

Note that s^Y^Li sm ( a j) > ^ £ n i s the area of the rhombic faces incident to 
the vertex z. Thus we conclude that /? = 0(e n ). This means, that if we choose 
j3 = Ce n with C > a sufficiently large constant and if e„ is small enough such 
that Ce n < 1, then p will have no maximum in Vi n t{G n )- In that case, as we 
have p(z) = {3\z\ 2 = 0(s n ) on Vo(G n ), we deduce that p(z) < 0(e n ) in F(G„) 
and thus 

t n (z)-h n {z) <0(e n ) for z€V(G n ). (28) 

The proof for the reverse inequality is almost the same. The only modifi- 
cations needed are reversing the sign of /? and a few inequalities. □ 

Remark 3.4. The statement of Lemma T3. 31 can be improved to 

h n (z)-t n (z)=0(e*) (29) 

in the case of a 'very regular' isoradial circle pattern. These are isoradial circle 
patterns such that for every oriented edge ej 1 = Zj x — z € E(G) incident to an 
interior vertex z € Vint (G) there is another parallel edge e j2 = Zj 2 — z € E{G) 
with opposite direction incident to z, that is e j2 = — e^. Furthermore, the 
corresponding intersection angles agree: a([z, ZjJ) = a([z, Zj 2 ]). This additional 
regularity property holds for example for an orthogonal circle pattern with the 



combinatorics of a part of the square grid, see Figure 2(a) 

The proof of estimation (j2"9"]) follows the same reasonings as above, but makes 
an C(e^)-analysis. The additional regularity implies that all terms of order 
vanish. 

Definition 3.5. For a function r\ : V(G) — * R define a discrete Laplacian by 
A V (z)= £ 2/i ([ , )aj]) (0)(T7(^) "»?(«))■ (30) 

[z,Zj]£E(G) 
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([z 1 ,z 2 ])(°) > one immediately has the following 

Lemma 3.6 (Maximum Principle). If Ar/ > on Vi n t{G) then the maximum 
of rj is attained at the boundary Vg (G) . 



The proof of Lemma 13.31 actually shows, that t n — h n is almost harmonic. 
More precisely, we have A(t n — h n ) = 0{e^). Adding a suitable subharmonic 
function (3\z\ 2 with (3 > big enough, we deduce that the resulting function p is 
subharmonic, that is Ap > 0, such that p attains its maximum at the boundary. 
This is also important for our proof of the following lemma. 



Lemma 3.7. Let t n and h n be defined as in Lemma Iff. 31 Let K C D be a 

compact subset in D. Then the following estimation holds for every interior 
vertex z € Vi n t{G n ) (~1 K such that all its incident vertices Z\, . . . ,Zi are also in 
V int (G n )nK: 

t n (zj) - K{ Zj ) - (t n (z) - h n {z)) = 0(e n (- log£ n )"5). (31) 

for j = 1, . . . , I. The constant in the O -notation may depend on K , but not on 
n or z. 

The proof of Lemma 13.71 uses the following estimation for superharmonic 
functions, which is a version of Corollary 3.1 of [SQ97] ; see also [SC971 Re- 
mark 3.2 and Lemma 2.1]. 

Proposition 3.8 ( SC97]). Let G be an undirected connected graph without 
loops and let c : E{G) — * R + be a positive weight function on the edges. Denote 
c(e) = c(x,y) for an edge e = [x,y] € E{G) and assume that 

sr^ c(x,z) 
m = max > — < oo. 

[x,y]£E!G) ' C(X,y) 

Denote by d(x,y) the combinatorial distance between two vertices x,y S V(G) 
in the graph G. Let B x (g) = {y G V{G) : d{x,y) < g} be the combinatorial ball 
of radius g > around the vertex x S V(G). Fix x € V{G) and R > 4 and set 

A R = sup g~ 2 W x (g), where W x {g) = ^ c(z,y). 

l - e - R zeB^ie), yeV(G) 

d(x,z)<d(x,y) 

Let u be a positive superharmonic function in B X (R +1), that is 

c(z, w)(u(z) — u(w)) < 

[z,w]GE(G) 

for all w £ B X {R +1). Let y be incident to x in G. Then 



u{x) 



u(y) 



Awi 2 JAr 

< = 

" vcjx, y) log 2 R 



Proof of Lemma \3. 7[ Lemma 13.31 implies that t n (zj) — t n (z) = 0(e n ) for all 
incident vertices z, Zj G V(G n ) since h n = log \g'\ is a C°°-function. Consider a 
Taylor expansion about of 

'J a (z^)(tn(Zj) ~ t n (z)) -TT. 
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Similar reasonings as in the proof of Lemma 13.31 imply that At n (z) = O(e^). 
Let p = t n — h n + (3\z\ 2 with (3 £ (0, 1). Choosing (3 = Ce n with a sufficiently 
large constant C > and e n small enough we deduce similarly as in the proof of 
Lemma EPl that Ap(z) > for all interior vertices z £ Vi n t(G n ). Now define the 
positive function p = e n + \\p\\ — p. Then Ap(z) < for all z £ Vi nt (G n ). The 
proof of Lemma 13. 31 shows that there is a constant Ci, depending only on g, D, 
and the labelling a, such that ||p|| < C\s n . Thus ||p|| < G?,E n with C 2 = 2C\ +1. 

To finish to proof, we apply Proposition l3.8l to the supcrharmonic function p. 
Remember that G n is a connected graph without loops and c(e) := 2/^, JO) > 
defines a positive weight function on the edges. The bound (TT5|) on the labelling 
a implies that 

c(x,z) 2tt cot(7r/4- C/2) 
m= max > — -< — — — — - — — — — =: 6 3 < oo. 

[x lV ]eE(G n ) [x z] f^ (Gn) c(x, y) tt/2 - C cot(7r/4 + C/2) 

Let x £ Vi n t(G n ). Note that 

W x (g) = V c(z,y)<( max c(e) ) \F w (x, g)\, 

d(x,z)<d(x,y) 

where F w (x, g) is the set of all faces of with one white vertex z £ B x (g) and 
\F w (x,g)\ denotes the number of faces of F w (x,g). Now, max e6 E(G„) c ( e ) < 
cot(7r/4-C/2)/2 and 

e ) C D x {{g + l)2e„) = {w £ C : |w - x\ < (g + l)2e„}, 

as the edge lengths in G n are smaller than 2e n . Remember that = 
e 2 sin «(/) > e 2 sin(7r/2 — C) is the area of the face / £ F{2) n ). Thus 



4 



„ 7T((g+ l)2e») 2 167rg : 
e^sm(7r/2 — C) sm(7r/2-6) 

for all g > 1. Therefore we obtain Ar = sup 1<e<R g~ 2 W x {g) < C4, where the 
upper bound C4 is independent of R > 4 and n £ IN. 

Let K C D be compact. Denote by e the Euclidean distance (between a 
point and a compact set or between closed sets of R 2 = C). Let z £ V(G n ) n X 
and set (i? + 1) = d(z, Vo(G n )) to be the combinatorial distance from z to 
the boundary of G n . Let Zj £ V(G„) be incident to z. As the labelling a is 
bounded, e„ — ► 0, and @ n approximates D, we deduce that R > e^Cg > 4 if 
n > no is large enough. Thus for all n > hq and all z £ V(G n ) n K 

1/v^og^R < l/V lo S2 - log 2 e n < V2/^-log 2 e n 
holds by our assumptions. Proposition 13.81 implies that 



^-c^z,,) log 2 ^ 



for all incident vertices z, Zj £ V(G n ) H if and n > no- As c(z, Zj) < cot(7r/4 
C/2)/2 and ||p|| < C 2 £n we finally arrive at the desired estimation 

\t n (zj) - h n (zj) - (t n (z) - h n (z))\ = \p(z) -p(zj)\ < C 6 e„(-log 2 £„)"5 
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for all incident vertices z, Zj G V(G n ) fl K and n > rig, where the constant Cq 
depends on C 2 , . . . , C5, that is only on g, D, C, C, and K. □ 

Lemma 3.9. Let e — uv G E(3l n ) be a directed edge with u G V{G n ) and 
v G V(G* n ). Denote by 5 n (e) the combinatorial distance in & n from e = [it, v] 
to [zq,vq], that is the least integer k such that there is a sequence of edges 
{[zq,vo] = ei, e2, . . . , ejt = e} C E{& n ) such that the edges e m+ i and e m are 
incident to the same face in H0 n for m = 1, . . . ,k — 1. Then 

tp n (e) = arg.g» + <f> n (e) + 5 n (e)0(e n (- log 2 e n )~*). (32) 

The constant in the notation 0(e n (— log 2 e n )~^) may depend on the distance 
of v to the boundary dD. 

Note that if 3D is smooth, then S n (e) = (^(e" 1 ). In general we have S n (e) = 
0(e~ 1 ) on compact subsets K C D, where the constant in the notation 0{e~ 1 ) 
may depend on K . In any case, on compact subsets of D we have 

y> n (g) = arg</ {v) + <j> n (g) + 0((- log 3 (33) 

Proof. Using the notation of Figure [3] (left), equation ([31) implies 

/ a (logr„(z+) - logr„(z_)) = f a (t n (z + ) - t n (z-)) 

= fa (0) + f a (0) (log \9'(z+)\ - log \g'(z-)\) + 0(e n (- log 2 e n )~i). 

As in Lemma l3~3l we have 2f' a (0) = |^~^~| = ^ with the same notation 

which yields 

2/;(0)(log| 5 '(z + )|-log| 3 '(z_)|) 

= .J + ~ t '', Re(a(z+ - *_)) + 0(e„(-log 2 e n )-') 
i(z+ - z_) 

= Im(a(v+ - + 0(e„(- log 2 £„)"') 

= ar g.9'(«+) - arg.g'(u_) + 0(e n (-log 2 e n )" 3 ), 
where a - s" ((z++z-)/2) _ fl "((„ + +„_)/2) 



By Lemma 12.71 we can choose the angle functions </> n and tp n on any mini- 
mal sequence of edges {[zq, vq] = ei, e 2 , . . . , = e} C E(@ n ) such that equa- 
tions ([T1)- ([T2"1) are satisfied without the (mod 27r)-term. Using the above consid- 
erations of 2/ Q (logr„(z + ) — log r n (z_)) and the normalization of ip n , we arrive 
at equation ([3"2"[) . □ 



Proof of Theorem \3.1\ Consider a compact subset K of D. Let z G V(G„) n if 
and v G U(G* ) n K be vertices which are incident in &! n , that is [z, v] G E(@ n ). 
Then Lemmas 13.31 and 13.91 imply that 



logs' (2) = log Iff' (2)| +iargg'(z) 

= \og(r n (z)/R n (z)) + i((p n (zv) - 4> n {zv)) + C((-log 2 £„)"'). 
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As g' and thus the quotient r r J R n is uniformly bounded, we obtain 

9 ' {Z) = ^tte*^-^"*"^" +©((- log 2 e n )-5). (34) 

S v ' 

=<!»(*) 

This implies the uniform convergence on compact subsets of D of g„ to g' . 

Convergence of g n is now proven by using suitable integrations of g' and q n . 
Let w € y(G* ) and consider a shortest path 7 in G* from «o to w with vertices 
{v = w u w 2 , ■ ■ ■ , w k = w} C F(G*). Then 

. fe-i 
g(w)=g(v ) + / g'{()d( = g(v ) + ^2g\w j+1 )(w J+1 - Wj ) + 0(e n ) 
J i 3=1 

k-1 

= 9M + 1n{wj+i)(w j+1 - Wj) + £>((- log 2 e n )-2 ), 
j'=i 

because g'(wj)-q n (wj) = 0((-log 2 e n )~5), w j+1 -Wj = 0(e n ) and k = ©(e^ 1 ) 
on compact sets. Thus it only remains to show that 

fe-i 

Pn(w) = g(v ) + Y 1n(w J + 1 )(w j+1 - Wj) + O ( (- log 2 £ n )~ 2 ) . (35) 
3=1 

Remembering 

«&.(«+) = ^4 f± T e ^" (SJTr ^ ) " " (, ^ TT ^ )) +0(("l0g a 

(w j+1 - w s ) = 2i?„(z + )cos(a(K- +1 ,^0)/2)c^"( s ^)-( 7r / 2 - a (^+ 1 '^])/ 2 », 

we can conclude that 

q n (Wj+i)(wj+i - Wj) =p n (wj+i) -p„(wj) + £>(£„(- log 2 e„)~^), 

where Z-,z+ € V{G) are incident to Wj+i and iu,' and we have used the notations 
are as in Figure [3] (left) with Wj = u_ and Wj + \ = v + . As we have normalized 
g( v a) = p(^o)) this proves equation (|35p and therefore the uniform convergence 
of p„ to g on compact subsets of D. □ 



Remark 3.10. Theorem 13.11 may easily be generalized in the following ways. 
First, we may consider 'nearly isoradial' circle patterns which satisfy R n (z) = 
0(e n ) for all vertices z € V{G n ) and i?„(zi)/i?„(z 2 ) = 1 + 0{e\) for all edges 
[z u z 2 ] €E(G n ). 

Second, we may omit the assumption that the whole domain D is approx- 
imated by the rhombic embeddings Sl n . Then the convergence claims remain 
true for compact subsets of any open domain D' C D which is covered or ap- 
proximated by the rhombic embeddings and contains vq. 

Using the angle function instead of the radius function, we obtain the fol- 
lowing analog of Theorem 13.11 for Neumann boundary conditions. 
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Theorem 3.11. Under the same assumptions as in Theorem \3.1\ and with the 
same notation, assume further that e n is sufficiently small such that for all 

sup max I &rgg'(v + 2e n e t0 ) — &xgg'{v)\ < — — C < min Ok — a{e)). (36) 

„ e .oee[o,27r] 2 e£E(G n ) 

Define an angle function on the oriented boundary edges by 

<Pn(e) = <f> n (e) + &rgg'(v), 

where e — zv £ E(2> n ) and v £ V{G* n ). Then there is a circle pattern ^ n for 
G n and a n with radius function r n and angle function ip n with these boundary 
values. 

Suppose that this circle pattern is normalized such that 

r n (z Q ) = R n (z )\g'(z )\, 

where zq £ V(G„) is chosen such that the disk bounded by the circle C ZQ contains 
0. Suppose further that ^ n is normalized by a translation such that 

PnM = g(v ), (37) 

where p n (v) denotes the intersection point corresponding to v £ V(G* ). 

Then q n — > g' and g n — > g uniformly on compact subsets in D as n — > oo. 

Proof. The existence claim for the circle pattern with Neumann boundary con- 
ditions follows from [BS04, Theorem 3] using the assumption (|3(il) . 

Theorem 12.81 shows that the difference ip n ~ <p n gives rise to a function 8 n : 
V{G* n ) — > R with boundary values given by argg'. The proof of the convergence 
claim is similar to the proof of Theorem 13.11 The roles of 5 n — ip n — <j) n and 
log(r„/i?„) have to be interchanged in Lemmas 13. 3[ 13. 7\ and 13.91 and similarly 
&rgg' = Im(logg') has to be considered instead of log \g'\ = Re(logp'). The role 
of equation ^ is substituted by equation (fT7|) . □ 



4 Quasicrystallic circle patterns 

The order of convergence in Theorems 13 . 1 1 and 13 . 1 1 1 can be improved for a special 
class of isoradial circle patterns with a uniformly bounded number of different 
edge directions and a local deformation property. In the following, we introduce 
suitable terminology and some useful results. 

As the kites of an isoradial circle pattern are in fact rhombi, an embedded 
isoradial circle pattern leads to a rhombic embedding in C of the corresponding 
b-quad-graph 'S . Conversely, adding circles with centers in the white vertices of 
a rhombic embedding and radius equal to the edge length results in an embedded 
isoradial circle pattern. 

Given a rhombic embedding of a b-quad-graph consider for each directed 
edge e £ E(@) the vector of its embedding as a complex number with length one. 
Half of the number of different values of these directions is called the dimension 
d of the rhombic embedding. If d is finite, the rhombic embedding is called 
quasicrystallic. A circle pattern for a b-quad-graph @ is called a quasicrystallic 
circle pattern if there exists a quasicrystallic rhombic embedding of & and if the 
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intersection angles are taken from this rhombic embedding. The comparison 
function of the isoradial circle pattern for S and the quasicrystallic circle 
pattern % will also be called comparison function for c &2- 

Quasicrystallic circle patterns were introduced in [BMS05] . Certainly, this 
property only makes sense for infinite graphs or infinite sequences of graphs 
with growing number of vertices and edges. 

In the following we will identify the b-quad-graph S> with a rhombic embed- 
ding of @>. 



4.1 Quasicrystallic rhombic embeddings and Z° 



Any rhombic embedding of a connected 
b-quad-graph S can be seen as a sort of pro- 
jection of a certain two-dimensional subcom- 
plex (combinatorial surface) f2@ of the multi- 
_/ dimensional lattice Z d (or of a multi-dimen- 
sional lattice C which is isomorphic to Z d ). 
' An illustrating example is given in Figure [5l 
The combinatorial surface in Z d can 
be constructed in the following way. Denote 
the set of the different edge directions of Si 
by A — {±ai, . . . ,±ad} C S 1 . We suppose 

that d > 1 and that any two non-opposite el- 
Figure 5: An example of a c . ,. , • j j, m 

, cmcnts ot A arc linearly mdcpcndent over JR. 
combinatorial surface S I <a C Z . T , , , , , 

Let ei, . . . , &d denote the standard ortnonor- 

mal basis of R d . Fix a white vertex xq & V{&) and the origin of R d . Add 
the edges of {±ei, . . . , ie^} at the origin which correspond to the edges of 
{±ai, . . . , ±cid} incident to xq in S, together with their endpoints. Successively 
continue the construction at the new endpoints. Also, add two-dimensional 
facets (faces) of Z d corresponding to faces of S, spanned by incident edges. 

A combinatorial surface f2@ in X d corresponding to a quasicrystallic rhom- 
bic embedding can be characterized using the following monotonicity property, 
see |BMS05( Section 6] for a proof. 

Lemma 4.1 (Monotonicity criterium). Any two points of VLo$ can be connected 
by a path in with all directed edges lying in one d- dimensional octant, 
that is all directed edges of this path are elements of one of the 2 d subsets of 
{±ei, . . . jie^} containing d linearly independent vectors. 

An important class of examples of rhombic embeddings of b-quadgraphs can 
be constructed using ideas of the grid projection method for quasipcriodic tilings 
of the plane; see for example [DK85 , GR86, Sen95j . 

Example 4.2 (Quasicrystallic rhombic embedding obtained from a plane). Let 
£ be a two-dimensional plane in R d and t £ E. Let ei, . . . ,ea be the standard 
orthonormal basis of JR d . We assume that E does not contain any of the seg- 
ments Sj = {t + \ej : A € [0, 1]} for j = 1, . . . , d. Then we can choose positive 
numbers ci, . . . ,Cd such that the orthogonal projections PE{cjGj) have length 
1. (If E contains two different segments Sj 1 and Sj 2 , the following construc- 
tion only leads to the standard square grid pattern Z 2 . If E contains exactly 
one such segment Sj, then the construction may be adapted for the remaining 
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(a) Straight line E. (b) Modified line E/\ with a "hat". 



Figure 6: Example for the usage of the grid- projection method in 1? . The 
Voronoi cells which contain Q c arc marked in red. 

dimensions excluding e,-.) We further assume that the orthogonal projections 
onto E of the two-dimensional facets Ej 1 ^ 2 = {Aie^ + \2^j 2 '■ Ai, A2 G [0, 1]} 
for 1 < ji < J2 < d are non-degenerate parallelograms. 

Consider around each vertex p of the lattice C = C\E X • • • X c<{L the 
hypercuboid V = [— ci/2,ci/2] x ••• x [— c<j/2, Qi/2], that is the Voronoi cell 
p + V. These translations of V cover R d . We build an infinite monotone two- 
dimensional surface SI (E) in C by the following construction. The basic idea 
is illustrated in Figure [H] (left). 

If E intersects the interior of the Voronoi cell of a lattice point (i.e. (p + 
V)° n E ^ for p G £), then this point belongs to Q C (E). Undirected edges 
correspond to intersections of E with the interior of a (d — l)-dimensional facet 
bounding two Voronoi cells. Thus we get a connected graph in C. An inter- 
section of E with the interior of a translated (d — 2)-dimensional facet of V 
corresponds to a rectangular two-dimensional face of the lattice. By construc- 
tion, the orthogonal projection of this graph onto E results in a planar connected 
graph whose faces are all of even degree (= number of incident edges or of in- 
cident vertices). A face of degree bigger than 4 corresponds to an intersection 
of E with the translation of a (d — A:)-dimensional facet of V for some k > 3. 
Consider the vertices and edges of such a face and the corresponding points 
and edges in the lattice C. These points lie on a combinatorial fc-dimensional 
hypercuboid contained in C. By construction, it is easy to see that there are 
two points of the fc-dimensional hypercuboid which are each incident to k of 
the given vertices. We choose a point with least distance from E and add it to 
the surface. Adding edges to neighboring vertices splits the face of degree 2fc 
into fc faces of degree 4. Thus we obtain an infinite monotone two-dimensional 
combinatorial surface il^(E) which projects to an infinite rhombic embedding 
covering the whole plane E. 

Example 4.3 (Modification of the construction in Example 14. 2p . The method 
used in the preceding example can be modified to result in similar but different 
rhombic embeddings. The basic idea is illustrated in Figure [S] (right). 

Let E be a two-dimensional plane in R d satisfying the same assumptions as 
in Examplc l4.2l Let N ^ be a vector orthogonal to E. Let A be an equilateral 
triangle in E with vertices ti,t2,t3 and let s be the intersection point of the 
bisecting lines of the angles. Consider the two-dimensional facets of the three- 
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dimensional tetrahedron T(A, N) spanned by the four vertices ti, t3, s + N. 
Exactly one of these facets is completely contained in E (this is the triangle 
A). We remove the triangle from E and add instead the remaining facets of 
T(A,N). Let -Ea be the resulting two-dimensional surface. Note that E/\ is 
orientable like E. Define 7 € (0,tt/2) by 7 = arctan(2v / 3||N||/||ti - t 2 |j), where 
|| • || denotes the Euclidean norm of vectors in III d . 7 is the acute angle between 
E and the two-dimensional facets of T(A, N) not contained in E. 

If 7 is small enough, then our assumptions imply that we can apply the same 
construction algorithm to E& as for the plane E in the previous example and 
obtain a monotone surface Q c (E/\). The orthogonal projection of Q c (E/\) onto 
E is a rhombic embedding which coincides with the rhombic embedding from 
il c (E) except for a finite part. 



The quasicrystallic rhombic embedding of Example 14.31 may also obtained 
using the following general concept of local changes of rhombic embeddings. 

Definition 4.4. Let ^ be a rhombic embedding of a finite simply connected 
b-quad-graph with corresponding combinatorial surface f2g> in r E d . 

Let z S Vint(Q@) be an interior vertex 
with exactly three incident two-dimensional 
facets of f2@. Consider the three-dimensio- 
nal cube with these boundary facets. Re- 
place the three given facets with the three 
other two-dimensional facets of this cube. 
This procedure is called a flip; see Figure [7] 
for an illustration. 

A vertex z £ I/ 1 can be reached with flips 
from fig, if z is contained in a combinatorial 
surface obtained from f2® by a suitable se- 
quence of flips. The set of all vertices which can be reached with flips (including 
V(Q@)) will be denoted by .F(f2®). 



Figure 7: A flip of a three- 
dimensional cube. The dashed 
edges and their incident faces are 
not part of the surface in 1i d . 



For further use we enlarge the set J- (fig,) of vertices which can be reached 
by flips from in the following way. For a set of vertices W C V(1i d ) denote 
by the set W together with all vertices incident to a two-dimensional facet 
of Z d where three of its four vertices belong to W. Define H/l fc+1 l = (I^M)! 1 ! 
inductively for all k G IN. In particular, we denote for some arbitrary, but fixed 
k e IN 

^„(fi») = (^(n»)) w . 



4.2 Quasicrystallic circle patterns and integrability 

Let *3> be a quasicrystallic rhombic embedding of a b-quad-graph. The combi- 
natorial surface in ifi is important by its connection with integrability. Sec 
also |BS08j for a more detailed presentation and a deepened study of integra- 
bility and consistency. 

In particular, a function defined on the vertices of which satisfies some 
3D-consistent equation on all faces of 11® can uniquely be extended to the brick 



:= {n = (m, ...,n d )eZ d : a k (tt®) < n k < 6 fe (fi@), k = l,...,d}, 



22 



Ulrike Bucking 



where a k (tt@) = tain ne v(o 9 ) n k and b k (fls>) = max ngV (nsi) n fc . Note that 
IT(f2@) is the hull of f2@. A proof may be found in [BMS05, Section 6]. This ex- 
tension of a function using a 3D-consistent equation will now be applied for the 
comparison function w defined in (fl4"|) of two circle patterns. In particular, we 
take for the isoradial circle pattern which corresponds to the quasicrystallic 
rhombic embedding & ', Given another circle pattern ^ for Q> with the same 
intersection angles, let w be the comparison function for ^2- Note that the Hi- 
rota Equation (fT5| is 3D-consistent; see Sections 10 and 11 of |BMS05j for more 
details. Thus w considered as a function on V(Q,&) can uniquely be extended 
to the brick II(f2@) such that equation (fT5|) holds on all two-dimensional facets. 
Additionally, w and its extension are real valued on white points V W (Q,@) and 
has value in ffi 1 for black points V&(f2@). This can easily be deduced from the 
Hirota Equation (fT5j) . 

The extension of w can be used to define a radius function for any rhombic 
embedding with the same boundary faces as S>. 

Lemma 4.5. Let @ and S>' be two simply connected finite rhombic embeddings 
of b- quad- graphs with the same edge directions. Assume that & and & agree on 
all boundary faces. Let be an (embedded) planar circle pattern for 3$ and the 
labelling given by the rhombic embedding. Then there is an (embedded) planar 
circle pattern c € l for <£l' which agrees with c € for all boundary circles. 

Proof. Consider the monotone combinatorial surfaces fl® and f2^. Without loss 
of generality, we can assume that Q@ and have the same boundary faces in 
X d . Thus they both define the same brick II(f2@) = Il(fi^) =: II. Given the cir- 
cle pattern ff, define the function w on V(£l@) by (fT4]l . Extend w to the brick II 
such that condition (fT5|) holds for all two-dimensional facets. Consider w on ClL 
and build the corresponding pattern c €' , such that the points on the boundary 
agree with those of the given circle pattern c 1o . Equation ijTaj) guarantees that 
all rhombi of fi^ are mapped to closed kites. Due to the combinatorics, the 
chain of kites is closed around each vertex. Since the boundary kites of c €' are 
given by which is an immersed circle pattern, at every interior white point 
the angles of the kites sum up to 2tt. Thus is an immersed circle pattern. 

Furthermore, c €' is embedded if is, because is an immersed circle 
pattern and ^f' and ^ agree for all boundary kites. □ 

5 C°°-convergence for quasicrystallic circle pat- 
terns 

In order to improve the order of convergence in Theorem 13.11 we study partial 
derivatives of the extended radius function using the intcgrability of the Hirota 
equation (|15p and a Regularity Lemma 15.61 

The following constants are useful to estimate the possible orders of partial 
derivatives for a function defined on JF«(Q®). Note that T K (yigi) C n(f2@). 

Definition 5.1. Let & be a rhombic embedding of a finite simply connected 
b-quad-graph with corresponding combinatorial surface Q@ in Z d . Let J C 
{1, . . . , d} contain at least two different indices. 
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For B > define a combinatorial ball of radius B about z G V(1i d ) using 
the directions {e.,- : j 6 J} by 

Uj(z, B) = {( = z + Y / n jGj : £ |n,-| < B}. (38) 

The radius of the largest ball about z using these directions which is contained 
in T K (n & ) is denoted by Bj(z, = max{B G IN : Uj(z,B) C ^(ftg,)}. 

Denote by d(z, dVL@) the combinatorial distance of z € V^(f2@) to the bound- 
ary <9f2@, that is the smallest integer K such that there is a connected path with 
K edges contained in from z to a boundary vertex of 9f2@. For further use, 
we define the constant 

Cj{F K <to 3 )) = min ( gj(Z '^io g) , ) + 1 : * G Vi„ t (fi») j > 0. (39) 
L d(z,di2@) J 

Note as an immediate consequence that for all z G Vi„{(0@) 

c/j(z, rc J (^ ( (n»))d(i I 5n») - 1]) c ^„(n»), 

where |~s] denotes the smallest integer bigger than s G R. 

The following theorem is an improved version of Theorem [3J] for a specified 
class of quasicrystallic circle patterns. 

Theorem 5.2. Under the assumptions of Theorem \3.1\ and with the same no- 
tation, let d G IN with d > 2 be a constant and assume further that S> n is a qua- 
sicrystallic rhombic embedding in D with edge lengths e n and dimension d n < d. 
The directions of the edges are elements of the set {±aj , . . . , ±a^™'} C S 1 such 

that any two of the vectors of {a\ , ...,a d > } are linearly independent. The 
possible angles are uniformly bounded, that is for all n G IN the scalar product 
is strictly bounded away from 1, 

|(a l (n) ,af ) )|<cos(V2 + C)<l, (40) 

for all 1 < i < j < d n and some constant < C < n/2. Consequently, the 
intersection angles a n are uniformly bounded in the sense that for all n G IN 
and all faces f G F{S> n ) there holds 

|a n (/)-7r/2|<C. (41) 

Let k G IN, let Jq C {1, . . . , d} contain at least two indices, and let B, Cj a > 
be real constants. Suppose that Jq C {1, . . . , d n } for all n G IN and 

Cj {F K {n® n ))>c Jo >o. 

Then we have with the same definitions of r n , <j) n , q n , and p n as in Theo- 
rem \3.1\ that q n — > g' and g n — > g in C°°(D) as n — > oo, that is discrete partial 
derivatives of all orders of q n and g n converge uniformly on compact subsets to 
their smooth counterparts. 

Simple examples of sequences of quasicrystallic circle patterns for this theo- 
rem are subgraphs of the suitably scaled infinite regular square grid or hexago- 
nal circle patterns or subgraphs of suitably scaled infinite rhombic embeddings 
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constructed in Examples 14.21 and 14.31 (see Figure \2§ ■ Simply connected parts 
of these rhombic embeddings which are large enough satisfy the conditions 
Cj{F K (£l@„)) > Co > for all subsets J c {l,...,d}, where the constant 
Co only depends on the construction parameters. This is a consequence of the 
simple combinatorics or of the modified construction in Example 14.31 

Remark 5.3. Similarly as for Theorem 13.11 the proof of Theorem 15.21 shows 
that we have in fact a priori bounds in 0(e n ) on compact subset of D for 
the difference of directional derivatives of log g' (and thus of g' and g) and 
corresponding discrete partial derivatives of log j 1 - + i((p n — <j) n ) (and thus q n 
and g n respectively). The constant in the estimation depends on the order of 
the partial derivatives and on the compact subset. 

Furthermore, for the 'very regular' case considered in Remark l3.4l the a priori 
bounds on the above partial derivatives have order 0(e„) on compact subset of 
D due to the improved estimation |2^|) . 

Remark 5.4. There is an analogous version of Theorem l5.2l of C°°-convcrgcncc 
for quasicrystallic circle patterns with Neumann boundary conditions. 

For the proof of Theorem 15. 21 we first define the comparison function w n for 
the circle pattern ^„ according to (|14[) and extend it to .F re (f2@ n ). This extension 
is again denoted by w n . The restriction of w n to white vertices is called extended 
radius function and denoted by r n as the original radius function for We 
also extend \g'\ to JF re (0® n ) by defining |</(z)| to be the value |</(z)| at the 
projection z of z <G Z d onto the plane of the rhombic embedding 3> n . If n is big 
enough, which will be assumed in the following, then z £ D or z £ W \ D and 
the distance of z to D is bounded independently of n. 

Denote t n = log(r n /i?„) = log(r„/e„) and h n = log|g'| as in Section [3] 
Using the extensions of r n and \g'\, these functions are defined on all white 
vertices of J- K (£l@ n ). Furthermore wc have the following extension of Lemma l3.3l 

Lemma 5.5. The estimation 

h n (z)-t n (z)=0(e n ). 

holds for all white vertices z £ V(J- K (£l@ n )). The constant in the O-notation 
may depend on k and on the dimension d n of S! n . 

Proof. Let z £ V{J-(£l$ n )) be a white vertex. By definition of there 
is a combinatorial surface f2'(z) containing z with the same boundary curve as 
VL^ n . Lemma 14.51 implies that we can define an embedded circle pattern %f' 
using the values of w n on f2^ n . Now the claim follows from Lemma 13.31 

For z £ V A (J r K (fJ S ) n ))\]/(J r (ri@ ii )) observe that equation (TT5|) may be used to 
extend the estimation of h n — 1„ . Each steps adds an error of order e n , therefore 
the final constant depends on k and on d n . □ 

Our main aim is to estimate the partial derivatives of the extended radius 
function in Z d " . Such partial derivatives can generally be considered in direction 
of the vectors v = ie^ ± ej 2 for < ji, j2 < d n such that and ej 2 are not 
collincar. Let vi, . . . , v 2 d n (d n -i) be an enumeration of these vectors and set 
V n = {vi, . . . , v 2( j, 1 (d„-i)}- The corresponding enumeration of the directions 
v = ±av^ ± aj T 2 l) in % for < j x ,j 2 < d is denoted by v x ,... V 2 d n {d n -i)- 



C°°-CONVERGENCE FOR QUASICRYSTALLIC CIRCLE PATTERNS 



25 



For any function h on white vertices z G Z d and/or in z G C define discrete 
partial derivatives in direction Vj or Vi by 

d Vi h(z) = : — : and o v ,h(z) = : : 

£n\Vi\ £n\Vi\ 

respectively. Furthermore, we call a direction v; or Vi to be contained in 
or 2t n at a vertex z or z respectively if there is a two-dimensional facet of fl @ n 
incident to z whose diagonal incident to z is parallel to v.;, that is {z + Av^ : 
A G [0,1]} C Q.@ n . 

Corresponding to these partial derivatives we use a scaled version of the 
Laplacian in (|30p . For a function 77 and an interior vertex z with incident 
vertices z±, . . . ,zl define the discrete Laplacian by 

A^n(z) := A-,..., V ^M :=iE 2 4(W)(°)(^i) " »?(*))■ ( 42 ) 

Here = Zj]) = (zj — z)/e n (j — 1, . . . , L) and the notation A^ v ^ 
emphasizes the dependence of the Laplacian on the directions u([z,2j]) of the 
edges [z, Zj] G E(G n ). 

Let zo G Vint(G n ) be an interior vertex and let v Ml , . . . , v Mi G V„ be the 
directions which correspond to the directions of the edges of G n incident to zq. 
Let zi G ^-" K (0@ n ) be a vertex such that Zi +v Mi G .^(f^) for all i = 1, . . . , L. 
Our next aim is to study A^™ v ^ t n (zi). For this purpose we assume that we 
have translated to zi the facets of f2^ n incident to Zo, that is we consider the 
(very small) monotone surface consisting of the two-dimensional facets incident 
to zi which contain a diagonal {zi + Av^ : A G [0, 1]} for i = 1, . . . ,L. Using 
the extension of the comparison function w n , the closed chain of these two- 
dimensional facets incident to Zi is mapped to a closed chain of kites. Thus we 
have 

L 

£ 2/a Mi (<n(*i + v Ml ) - tnfa)) G 27T1N, 
i=i 

where denotes the labelling of the two-dimensional facet containing Zi and 
zi + v w . By assumption X)i=i ^/a^ (0) = 27r and we know that t n (zi + v w ) 



t n (zi) = C(£ n ) by Lemma 15.51 Since the intersection angles and thus the 
maximum number of neighbors arc uniformly bounded, we deduce that 

{j2 2 ^ (^ Zl + V ^ ) - - 27T = 

if e n is small enough. Using a Taylor expansion about 0, we obtain 

9 ( L 00 f^fO) \ 
A- i ,..., v ^t„(z 1 ) = tj -EE ^(^i+v B )-UC 

™ \ 1=1 m=3 ' / 

/ L 00 /-(m) ("n") \ 

= e « h 2 ^ S ^Kr(^(zorc- 3 

\ i=l m=3 ' / 

— : £nF Vltl ,...,v Mi (^ni ^v M1 t n , . . . , d VfiL t n ; zi), (43) 
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Note that F v v is a C°°-function in the variables e„, d v t n , . . . , d v t n . 
This fact will be important for the proof of Lemma [5771 below. 

For further use we introduce the following notation. Let K be a compact 
subset of D and let M > 0. Denote by f2^' M the part of fl%> n with vertices of 
combinatorial distance bigger than M to the boundary and whose corresponding 
vertices z G V(@ n ) lie in K. Let J C {1, . . . ,d} contain at least two different 
indices. In order to consider partial derivatives within K in the directions 
iejj ± e j2 , where ji, ji G J and ji we attach a ball Uj(z, M) at each of 

these points: 

Uj(K,M,H 9n )= |J Uj(z,M). (44) 

Note that if M < Bj(z), for example if M < Cj(T K {Q.® n ))d(z, <9S1@„) - 1 for 
all z G V(n*f), then Uj(K,M,Q® n ) C ^(fiflj. 

Furthermore we define AT + do to be the compact do-neighborhood of the 
compact set K C C, that is 

K + do = {ze€:e(K,z)<do}, 

where <e denotes the Euclidean distance between a point and a compact set or 
between two compact subsets of C. 

The following lemma is important for our argumentation, as it gives an esti- 
mation of a partial derivative using estimations of the function and its Laplacian. 

Lemma 5.6 (Regularity Lemma). Let & be a quasicrystallic rhombic embedding 
with associated graph G and labelling a. Let W C V(G) and let u : W — > R be 
any function. Let M{u) = max^giy^j \ Au(v) / (4F* (v))\, where 

F*(v) = \ £ c{[zM)\z-v\ 2 = \ J2 

[z,v]£E(G) [z,v]eB(G) 

is the area of the face of the dual graph G* corresponding to the vertex v G 
Vi n t(G). There are constants C^^Cq > 0, independent ofW and u, such that 

\u(x ) - u( Xl )\p < C 5 \\u\\ w + p 2 C 6 M{u) (45) 

for all vertices X\ G W incident to Xq, where p is the Euclidean distance of Xq 
to the boundary Wq. 

A proof is given in the appendix, see Lemma IA.81 As a direct application 
we get 

Lemma 5.7. Let K C D be a compact set, no € IN and < do < e(K,d3t n ) 
for all n > no ■ 

Let k G Mo and let , . . . , Wi k G Hn>n ^ n ^ e ^ no ^ necessa7 "dy differ- 
ent directions. Let J C {1, ...,d} be a minimal subset of indices such that 
{v^, . . . , Vi fc } C span{ej : j G J}. Let Bq, Cq > be some constants. Assume 
that all discrete partial derivatives using at most k of the directions , . . . , Vi k 
exist on Uq = Uj(K + do, Bo,Q@ n ) and are bounded on Uq by C§e n for all 
n>n$. 

Let n > uq and let f2^ be a two-dimensional monotone combinatorial sur- 
face. Let 2> n be the corresponding rhombic embedding with edge lengths e n and 
such that 

{z g v(si § )-.ze v($ n ) n(K + d )} c u . 
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Let Zq £ V w (@ n ) fl (K -\-do/2) be a white vertex such that <e(z, dS> n ) > do jl. Let 
v *fc+i G V &e a direction contained in fig, at zq. Then there is a constant C\, 
depending on K , D, g, and on the constants do, Bo, Co, k, C , but not on zq, 
Vi k+1 , and n, such that 

\d Vtk+1 d Vik ■ ■■d Vii (t n - h„)(z )\ < C x e n . 

Proof. The proof is an application of the Regularity Lemma l5.6l for the function 
u = d Vik ■ ■ ■ d Vii (t n — h n ) and the part of the rhombic embedding 3> n contained 

in K + do- Note that the edge lengths of £P„ are e n and we suppose that 
*(z,d9 n ) > do/2. 

By assumption we have v(@ n )n(K+d ) — ^o^n- As h n is a C°°-function, 
this implies that d Vi d Vik ■ ■ ■ d Vii t n = 0(1) on Uq for all possible directions v, £ 
V whenever this partial derivative is defined in T K (^l^ ). 

Let Z\ £ Vi n t{G n ) n (K + d ) be an interior white vertex of 3> n and let 
v Ml , . . . , v ML £ V„ be the directions corresponding to the directions of the 
edges of G n incident to z\ . Note that 

Kl 1 ,..., v ^= Kl^..., v . L d Vl -- -d Vn (t n -h n ) = d Vt -- -d Vii A^ it ^ tV ^ (t n -h n ). 

From the above consideration in (1431) we use that A5" t n = e n F v v , 

where Fv M is a C°°-function in the variables e n , 9 V(li t n , . . . , 9 Vf> t n . From 

our assumptions we know that all partial derivatives <9 Vm . <9 Vifc ■ • • d Vii t n (zi) and 
those containing less than fe+1 of the derivatives (9 V(J . , <9 Vifc , ■ ■ ■ , d Vii arc defined 
and uniformly bounded by a constant independent of z\ and e n . Thus 

— Cv lil ,...,v l t L £n, 

where the constant C Vmi .„ iV(1 does not depend on z\ and e n . As h n is a 
harmonic C°°-function, we obtain by similar reasonings as in the proof of 
Lemma 13.31 that 

d Vik ---d Vii A^ it _ VhL h n = 0(e n ). 

As V„ is a finite set, we deduce that A e ™u is uniformly bounded on Vi„t(G„) fl 
[K + do) by C^En for some constant C2 independent of e n . Now the Regularity 
Lemma 15.61 gives the claim. □ 

The following corollary of the proceeding lemma constitutes the crucial step 
in our proof of Theorem 15.21 

Lemma 5.8. Let K C D be a compact set and let < d\ < e(K,dD). Let 
no £ W be such that K + d\ is covered by the rhombi of S# n for all n > no ■ Then 
there is a constant C\ = C\{K, C j^T^Q,®^))) > such that for all z £ K + d\ 

Cie- 1 < C Jo (F K (fl® n )) ■ d{z,d% x ) - 1. 

Furthermore, let k £ JNo and let , . . . , v; fc £ V be k ( not necessarily 
different) directions such that £ span{ej : j £ Jo} for all I — 1, . . . , k. Then 
there are constants no < n\(k,K) £ IN and C(k,K) > which may depend on 
k, K, d\, D, g, k, Cj (.^(f^J), but not on e n , such that for all n > n\(k,K) 
we have 

\\d Vik ---d Vii {t n -h n )\\ Uk <C{k,K)e ni (46) 
where U k = U Jo (K + 2-%, 2~ k C 1 s^\ fi^J. 



d ^ k '-- d ^ 1 A Z 1 ,...,v tl Jn(zi) 
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Proof. The existence of the constant C\ follows from the fact that d(z, d@ n )/£n 
is bounded from below for z £ (A' + d{) since the distance e(A + di,d@ n ) > 
is positive and the angles of the rhombi are uniformly bounded. 

The proof of estimation (jUJ) uses induction on the number of partial deriva- 
tives k. For k = the claim has been shown in Lemma 15.51 

Let k £ Mo and assume that the claim is true for all v < k. Let Vi fc+1 £ V 
be a direction with ~Vi k+1 £ {^ e ji i e j 2 } C spanje^- : j £ Jo}. Using the 
induction hypotheses, we can apply Lemma 15.71 for Uq = Uk, do = 2~ k d\, 
ileg — U{j 1 j 2 }(K + 2~ k di, 2 _fc Ci£~ 1 , f2g>„) C Uk and the corresponding rhom- 
bic embedding $l n obtained by projection, and z e!?„n (K + 2~ k ~ 1 d 1 ). This 
completes the induction step and the proof. □ 

Proof of Theorem \5.Sl Identify C with R 2 in the standard way and fix two or- 
thogonal unit vectors ei, e 2 - Define discrete partial derivatives d ei , d e2 in these 
directions using the discrete partial derivatives in two orthogonal directions 
Vi 1 = Oj^ + and Vi 2 = — a^ 1 for ji, .72 £ Jo- This definition depends 

on the choice of aj^ , a^™ , which may be different for each n, but this does not 
affect the proof. As the possible intersection angles arc bounded and as h n is a 
C°°-function, we deduce that 

\\d ejk ■ --dentin - d jk ■ ■ ■ d^hnWn < Ci(k,K)e n 

on every compact set A for jk, ■ ■ . , ji £ {1,2}. Here d\, d 2 denote the standard 
partial derivatives associated to ei,e2 for smooth functions and Ci(fc, A) is a 
constant which depends only on A, k, and <?. Lemma 15.81 implies that 

ll^ e i fc ■ ■ ■ 9e n t n - de jk ■ ■ ■ de^hnWu hi ( K+2 -k du 2-k Cie -^ < C 2 (k,K)e n 

if n is big enough. Using a version of Lemma 13.91 with error of order C(e„), we 
deduce that t n +i(tp„ — </>„) converges to \ogg' in C°°(D). Now the convergence 
of q n and g n follows by similar arguments as in the proof of Theorem 13.11 □ 

A Appendix: Properties of discrete Green's func- 
tion and regularity of solutions of discrete el- 
liptic equations 

In order to prove the Regularity Lcmma l5.6l (sec Lemma IA.8|) . we present some 
results in discrete potential theory on quasicrystallic rhombic embeddings which 
are derived from a suitable asymptotic expansion of a discrete Green's function. 

Throughout this appendix, we assume that S> is a (possibly infinite) simply 
connected quasicrystallic rhombic embedding of a b-quad-graph with edge di- 
rections A = {±ai, . . . ,±a,d}- Also, the edge lengths of & are supposed to be 
normalized to one. Let G be the associated graph built from white vertices. 

Fix some interior vertex xq £ Vi n t(G). Following Kenyon [Ken02| and 
Bobenko, Mercat, and Suris |BMS05| . we define the discrete Green's function 
g(x , •) : V(G) -> R by 

g(x , x) = - J^ jf l ^jp-e(x; X)d\ (47) 
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for all x G V(G). Here e(x; z) — Y\k=i y z-at j * s ^ nc discrete exponential 
function, where n = (m, ■ ■ ■ , ^d) = x — xo G Z d and x, Xo € 1^(0®) correspond 
to x,xq G V(S>) respectively. The integration path V is a collection of 2c? 
small loops, each one running counterclockwise around one of the points zta^ 
for k = 1, . . . , d. The branch of log(A) depends on x and is chosen as follows. 
Without loss of generality, we assume that the circular order of the points 
of A on the positively oriented unit circle S 1 is a\, . . . , ad, -o.it ■ ■ ■ , —ad- Set 
o-k+d = —a-k for k = 1, . . . , d and define a m for all m G Z by 2d-periodicity. To 
each a m = e m G S 1 we assign a certain value of the argument m G R: choose 
6\ arbitrarily and then use the rule 

Om+i - 8m G (0, 7r) for all m G Z. 

Clearly we then have Q m +d — 9 m -\-Tt. The points a m supplied with the arguments 
8 m can be considered as belonging to the Riemann surface of the logarithmic 
function (i.e. a branched covering of the complex A-plane). Since is a mono- 
tone surface, there is an m G {1, . . . , 2d} and a directed path from xq to x in 
@ such that the directed edges of this path are contained in {a m , . . . , a m +d-i}, 
see [BMS051 Lemma 18]. Now, the branch of log(A) in (|47p is chosen such that 

log(az) G [id m ,id m+ d-i], I = m, . . . , m + d - 1. 

Remember Definition 13. 51 of the discrete Laplacian and the representation of 
its weights c([zi,z 2 ]) = 2/^,, Z2 i)(0) in j4}. Then there holds 

Lemma A.l ( Kcn02, Theorems 7.1 and 7.3]). The discrete Green's function 
G(xo,-) defined in equation (|47[) has the following properties. 

(i) AQ(xo,v) = — S Xo (v), where the Laplacian is taken with respect to the 
second variable. 

(a) G(x ,x ) = o. 
(Hi) g(x ,v) = O(log(\v-x \)). 

Note that Q(xq, •) may also be defined by these three conditions. 

A.l Asymptotics for discrete Green's function 

Kenyon derived in |Ken02| an asymptotic development for the discrete Green's 
function using standard methods of complex analysis. His result can be slightly 
strengthened to an error of order 0(l/\v — xo\ 2 ). Note, that there is the sum- 
mand — log2/(27r) missing in Kenyon's formula (but not in his proof). 

Theorem A.2 (cf. |Ken021 Theorem 7.3]). For v G V(G) there holds 

SM - -I log,2|„ - 10 |) - 2ffl + O (,-Lp) . (48) 
Here ^Euier denotes the Euler 7 constant. 
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Proof. Consider a directed path xq = Wq, . . . ,Wk = v in S from xq to v such 
that the directed edges of this path are contained in {a rra , . . . , a rn +d-i\ for some 
1 < m < 2d as above. Note that k is even since xq and v are both white vertices 
of $). The integration path T in (|^7|) can be deformed into a connected contour 
lying on a single leaf of the Riemann surface of the logarithm, in particular 
to a simple closed curve Ti which surrounds the set {a m , . . . , a m +d-i} in a 
counterclockwise sense and has the origin and a ray = {se l9 : s > 0} in its 
exterior. We also assume that T\ is contained in the sector {z = re tv : r > 
0, (p G + f + ^ — 77]} for some 77 > independent of w, xo, and m. This 
is possible due to the fact that 9 m+ d-i — S m < n — S for some S > independent 
of v, ceo, and m. 

Let N = \v — Xo\. Take < Q\ >C and g 2 ^$> N 3 , but not exponentially 

smaller than 1/iV or bigger than N respectively. The curve T% is again homo- 
topic to a curve T2 which runs counterclockwise around the circle of radius Q2 
about the origin from the angle to 8 + 2ir, then along the ray Si from Q2B 
to Qie , then clockwise around the circle of radius q\ about the origin from the 
angle 8 + 2tt to 8, and finally back along the ray from gie l9 to ■ Without 
loss of generality, we assume that M is the negative real axis. 

Kenyon showed in [Ken02 that the integrals along the circles of radius Q\ 
and Q2 give 

(-1)^(1 + 0(N Q1 )) - ^(1 + 0(N/ Q2 )). 

The difference between the value of log z above and below the negative real 
axis is 2iri. Thus the integrals along the negative real axis can be combined into 



11 1 z + bj 
7 II 



82 j = Q 3 



where bj — ifj+i — wj € A is the directed edge from Wj to Wj+i and fc = O(N) 
is the number of edges of the path. This integral can be split into three parts: 
from — Q2 to — y/N, from —\f~N to —l/y/N, and from — 1/yN to —Q\. 
The integral is neglectible for the intermediate range because 



t + e i/3 



< e 



It cos/3/(t-l) 2 



for negative t < and due to our assumptions. 
For small Itl we have 



1 t + ^ t-r tbj + 1 



II j-r. - II ir^r = (-^^ ^(i + oc** 3 )), 

i=o 3 3=0 1 

using the Neumann series and a Taylor expansion. Thus the integral near the 
origin is 

(— \\ k ( r~ ei p2(v-x )t r-Qi 2(v-x )t \ 

_LiL / ! dt+ o(kt 3 )- dt). 

4?r \J-i/Vn 1 J-i/Vn t J 
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Applying similar reasonings and estimations as in Kenyon's proof, we obtain 
(-l) fc ( 1 

(l0g(2£l(« - Xq)) + 7Eulcr) + O ' 



Here 7Euier denotes the Euler 7 constant. 

For large \t\ the estimations are very similar. Since 

we get 

log f - g2 + 7 Eulcr) + O ' 



Since fc is even, the sum of all the above integral parts is therefore given by 
the right hand side of ([48]). □ 



For a bounded domain we also define a discrete Green's function with vanish- 
ing boundary values. Let W C Vi n t(G) be a finite subset of vertices. Denote by 
Wg C W the set of boundary vertices which are incident to at least one vertex 
in V(G) \ W. Set W mt =W\Wg the interior vertices of W. Let x Q € W mt be 
an interior vertex. The discrete Green's function Gw(xq, •) is uniquely defined 
by the following properties. 

(i) AGw{xo, v) = —S Xo (v) for all v € Wmt, where the Laplacian is taken with 
respect to the second variable. 

(ii) Gw{x , v) = for all v e W 9 . 

In the following, we choose W to be a special disk-like set. Let xq € V(G) 
be a vertex and let p > 2. Denote the closed disk with center xq and radius 
p by B p (xq) C C. Suppose that this disk is entirely covered by the rhombic 
embedding @. Denote by V(xq, p) C V(G) the set of white vertices lying within 
B p (xo)- For x\ € Vi n t(xo,p) we denote 

Gx ,p{xi, ') = Qv{x ,p){Xl, ■)■ 

The asymptotics of the discrete Green's function G from Theorem IA.2I can 
be used to derive the following estimations for Gx ,p- 

Proposition A. 3. There is a constant C\, independent of p and Xq, such that 

\Gx , P {x ,v)\ < Ci/p 

for all vertices v € Vint (^0 1 P) which are incident to a boundary vertex. 

Furthermore, there is a constant C2, independent of p and Xq, such that for 
all interior vertices x\ <E Vi„t{xQ,p) incident to xq and all v <E V(xo,p) there 
holds 

\Gx ,p{xo,v) - Gx a ,p(xi,v)\ < C2/O - X \ + 1). 
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Proof. Consider the function h p (xo, ■) : V(xq,p) — ► R defined by 

h p (x ,v) = G Xa ,p{x ,v) - G(x , v) - ^(l0g(2p) + 7Eulor)- 

Then h p (xQ, •) is harmonic on Vint (%o , p) ■ For boundary vertices v € Vd(xo,p) 
Theorem IA.2I implies that h p (xo,v) — 0(1/ p). The Maximum Principle 13.61 
yields \h p (xo, v)\ < C/p for all v G V(xo,p) and some constant C independent 
of p and v. This shows the first estimation. 

To prove the second claim, we also consider the harmonic function 

h p (xi,ll) = Q Xo ,p{x\,v) -G{X!,V) - ^-(l0g(2p) +7Eulor) 

for a fixed interior vertex x\ G Vi„t(xo, p) incident to xq. By similar reasonings 
as for h p (xo,-), we deduce that \h p (x\,v)\ < C/p for all v G V(xo,p) and 
some constant C independent of p and v. Theorem IA.2I implies the desired 
estimation. □ 



A. 2 Regularity of discrete solutions of elliptic equations 

In the following, we generalize and adapt some results of discrete potential 
theory for discrete harmonic functions obtained by Duffin in |Duf53[ see in 
particular Lemma 1 and Theorem 3-5] . The proofs are very similar or use ideas 
of the corresponding proofs in |Duf53| . Our aim is to obtain the Regularity 
Lemma IA.8I 

Lemma A. 4 (Green's Identity). Let W C V(G) be a finite subset of vertices. 
Let u, v : W — > R be two functions. Then 

^2 (v(x)Au(x) - u(x)Av(x)) = ^2 c (,[p,q\)( v (p) u (q) - u (p) v (q)), 
xeWi nt \p,q]eE d (w) 

(49) 

where E a {W) = {\p, q] G E(W) :p&W mt ,qe W a } . 

Corollary A. 5 (Representation of harmonic functions). Let u be a real valued 
harmonic function defined on V(xo,p). Then 

q&Vg(x Q ,p) 

where 

c(q)= c([p,q})G X0 , p (x ,p) = O(l/p). (50) 

and [p,q]£E(G) 

The estimation in (|50|) is a consequence of Proposition lA.3l and of the bound- 
cdncss of the weights c(e). 

Theorem A. 6. Let u : V(xo,p) — * R be a non-negative harmonic function. 
There is a constant C3 independent of p and u such that 

V£Vint(xO,P) 

where F*(v) = \J2[ Z v]eE(G) c ([ z ' v ])\ z ~ v \ 2 * s area °f ^ e f ace °f ^ e dual 
graph G* corresponding to the vertex v G Vi n t(G). 



< 



C 3 u(x ) 



(51) 
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Proof. Consider the function 

p(z) = g( Xo , z) + ^-(iog(2 P ) + 7Eulcr ) + v ~ x f~ p2 . 

2ir Amp 1 

Let z £ Vi n t(G) be an interior vertex. Consider the chain of faces fi, ■ ■ ■ , f m of 
SI which are incident to z in counterclockwise order. The enumeration of the 
faces fj and of the black vertices v\ , . . . , v m and the white vertices z±, . . . , z m 
incident to these faces can be chosen such that fj is incident to Wj-i, Vj, and zj 
for j = 1, ... ,771, where vq = v m . Furthermore, using this enumeration we have 

zj - z . Vj-vj-i \vj-Vj-i\ 

-r — \ l = -r — — i i — -r-( z i- z ) = - 1 ); 

|^-^| \Vj-Vj-i\ \Zj-A 

see FigureUJwith z_ = z, z + = Zj, t>_ = Uj-i, t>+ = «j. As — a;o | 2 — \z — xq\ 2 = 
— 2Re((z — Xo)(Zj — z)) + \zj — z\ 2 we deduce by very similar calculations as in 
the proof of Lemma RIT31 that 



- m 

Ap(z) = Ag(x ,v) + ^—^Y^ cjfazj]) (\ z j ~ x o\ 2 - \z - x \ 2 ) 

= -5 Xo (z) + F*(z)/(np 2 ). 
Let v be incident to a vertex of Vq(xq, p). Theorem I A . 2 1 implies that 

P(v) = -I log + \ V -f-f + O ( T *) = 0(1/ p 2 ). 

2tt p 4np z \\v — xq\ 



Thus there is a constant B\, independent of p and v, such that px(v) := p(v) + 
Bi/p 2 > and |pi(w)| < 2B\j p 2 for all vertices v £ Vi nt {xo,p) incident to a 
vertex of Vg(xo,p). Applying Green's Identity lA. 41 to p\ and the non- negative 
harmonic function u, wc obtain 

u(x ) \ F*(v)u(v)= (pi(x)Au(x) - u(x)Api(x)) 

irp z z — ' — ' 

vev int (xo,p) xev int (xo,p) 

= T] c([z,q]){p 1 (z)u(q)-u(z)p 1 (q)), 

S V ' 

[z,q]£E p >Q 

^-4?r V u(q) < 8lrBlB2 u(x ), 
p z p 

q£V a (x ,p) 



< 



Here E p = {[x,y] £ E(G) : x £ Vi n t{xo, p), y £ Vg(xo,p)} and we have used 
the estimation 

]T c([x,y})< ]T c([x,y])\x-y\ 2 =4F*(x) < Air 

[x,y]GE(G) [x,y]GE(G) 

for all fixed vertices x £ Vi„t(G). Furthermore Z~2 v ^ v f ,(x n .p ) u (y) — ^2pu(xo) for 
some constant B2 > as a consequence of Corollary I A. 5l 

For the reverse inequality, note that there is also a constant B3 independent 
of p and v such that pi(v) := p(v) — B 3 /p 2 < and |pi(v)| < 2B 3 /p 2 for 
all vertices v £ Vi„t(xo, p) incident to a vertex in Vq{xq,p). Combining both 
estimation proves the claim. □ 
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Theorem IA.6I can be interpreted as an analog to the Theorem of Gauss in 
potential theory. Furthermore, we can deduce a discrete version of Holder's 
Inequality for non-negative harmonic function. 

Theorem A. 7 (Holder's Inequality). Let u : V(xo,p) — > R be a non-negative 
harmonic function. There is a constant C4, independent of p and u, such that 

\u(xq) - u{xi)\ < C 4 u(x a )/p (52) 

for all vertices X\ £ V(xo,p) incident to xo. 

As a corollary of Holder's Inequality and of Proposition IA.3I we obtain the 
following result on the regularity of discrete solutions to elliptic equations. 

Lemma A. 8 (Regularity Lemma). Let W C V(G) and let u : W — > R be any 
function. Set M{u) = max ve w int | / (4J^* (v) ) | , where F*(v) is the area of 

the face dual to v as in Theorem \A.6l Define ||t?||w : = max{|7y(2;)| : z £ W}. 
There are constants C$,C§ > 0, independent of W and u, such that 

\u(x ) - u(xi)\p < C 5 \\u\\ w + P 2 C 6 M{u) (53) 

for all vertices X\ £ W incident to xq G Wint, where p is the Euclidean distance 
of xq to the boundary Wq . 

Proof. Let x\ £ W be a fixed vertex incident to Xq. 

First we suppose that p > 4. Consider the auxiliary function f(z) = 
M(u)\z — xq\ 2 . Since \x± — xq\ < 2, we obviously have 

\f(x ) - f( Xl )\ = M(u)\x! - x \ 2 < 4M(u). 

Let h : V(xo, p) — > R be the unique harmonic function with boundary values 
h(v) = u(v) + f(v) for v £ Vg(xo, p). Holder's Inequality (f5"2")l and the Maximum 
Principle 13.61 imply that 

\h(x ) - h{xi)\p < Bi||/i||v(* 0l p) < Bidklk + M{u)p 2 ) 

for some constant B\ independent of h, p, x$, x\. 
Next consider s = u + f — h on V(xq, p). Then 

As = Au + AF*M(u) > on V int (x ,p), 
s(v) = for v £ Vg(xo,p). 

The Maximum Principle 13.61 implies s < 0. Green's Identity lA.41 gives 
s(x )+ ^2 G Xo .p(x ,v)As(v) 

veVint&OiP) 

= (Gx„A x o,v)^s{v) - s(v)Ag x0tP (x ,v)) 

v£V in t(x ,p) 

= X] c (iP^ ( l\)(Sx a ,p(xo,p)s(q) - s(p)g xo . p (xo,q)) = 0, 
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where E p = {[p, q] G : p G V$ nt (a:o, p), 5 € Va(x , p)}. Analogously, 

s(xx)+ ^ ^ 0! p(xi,w)As(u) = 0. 

Using the estimation As(v) < 8F*(v)M(u) we deduce that 
\s(x ) - s(xi)\ < 

weVi„t(xo,p) 



Now Proposition I A. 31 implies that 



\s(xo) - s( Xl )\ < 8B 2 M(u) Yl | r TT 1 ^ 8B 2 M(u)B 3 p, 

vev int (x , P ) lv x ° l + 

where B 2 and B% are constants independent of s, p, xq, x\. 

Combining the above estimations for f,h, and s, we finally obtain 

\u(x ) - u(xi)\p < \s(x ) - s(x-i) - (f(xo) - f(xi)) + h(x ) - h(x{)\p 
< B 1 \\u\\ w + p 2 (A + B 1 + 8B 2 B 3 )M(u). 

This implies the claim for p > 4. For p < 4 inequality (|53p can be deduced from 
-4F*(x Q )M(u) < Au(x ) = Y c([x ,v]){u(v)-u(x Q ))<4:F*(x )M(u). 

[x ,v]eE(G) 

using F*(xq) < ir, J2[ Xo v]eE(G) c ([ x 0t v }) < 47r, and the uniform boundedness of 
the weights c(e). □ 
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